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We non-perturbatively calculate the scale dependence of the static axial current in the 
Schrodinger functional scheme by means of a recursive finite-size scaling technique, taking 
the continuum limit in each step. The bare current in the O(a) improved theory as well as in 
the original Wilson regularization is thus connected to the renormalization group invariant 
one. The latter may then be related to the current at the B-scale defined such that its 
matrix elements differ from the physical (QCD) ones by 0(1/M). At present, a (probably 
small) perturbative uncertainty enters in this step. As an application, we renormalize ex- 
isting unimproved data on _p^ are and extrapolate to the continuum limit. We also study an 
interesting function h(d/L, u) derived from the Schrodinger functional amplitude describing 
the propagation of a static quark-antiquark pair. 
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1 Introduction 



Since the decay constant governing the leptonic decay of a B-meson, is an essential 
element in the quantitative analysis of the unitarity triangle, many lattice QCD investigations 
have worked towards its determination. However, with its large mass, the b-quark still 
escapes a direct numerical treatment [J and effective theories have to be used to separate 
the large mass scale from the low-energy bound-state dynamics. (As an exception to this 
rule, it has recently been demonstrated that also finite-volume methods on lattices with a 
large number of points represent a possible alternative 

The most natural and theoretically appealing effective theory is the static approxima- 
tion 4 . It describes the large-mass limit of the theory and is the starting point for a 
1/m-expansion, the heavy quark effective theory. Yet the problems of this framework have 
been twofold in the past, (i) The renormalization properties of the static theory are different, 
i.e. the renormalization constant of the axial current in (^ at )o = Z^ at (fj,) V'd7o75' i /'b tat 
becomes scale (/x) dependent, thereby entailing an additional uncertainty, and (ii) Monte 
Carlo computations of the matrix element itself are difficult. For these reasons, after a sig- 
nificant initial effort [51|Hl[7||H l|^lllL)lllllll2j . the compulation of in the static approximation 
has received little attention in the recent past. 

In the present work, we solve (i) by computing the renormalization factor that relates 
the bare operator in lattice regularization to the renormalization group invariant (RGI) 
operator. Denoting its matrix element by one then has a relation 

FbV^ = C PS (M h /A m ) x + 0(1/M b ) (1.1) 

with a function Cps of the renormalization group invariant b-quark mass Mb in units of 
the A-parameter. It is scale independent, but in practice it is obtained perturbatively and 
an uncertainty due to perturbation theory remains, see Section [21 The important task of a 
lattice computation of the B-meson decay constant in the static approximation is to compute 

^RGI- 

Our strategy to arrive at ^f^Qj from the bare matrix element follows the one used by 
the ALPHA Collaboration for the renormalization group invariant quark mass ^3]. In this 
approach, an intermediate finite-volume renormalization scheme is used to follow the scale 
evolution non-perturbatively to high energies (O(lOOGeV)), where then perturbation theory 
can safely be used to connect to the renormalization group invariants. For a more detailed 
explanation of the overall strategy we refer to Ref. ^1] and for a preliminary report on our 
work to Refs. [ToUTo]. 

The present paper is organized as follows. In Section |2] we discuss our intermediate 
renormalization scheme, formulated in the Schrodinger functional (SF). It is based on the 
renormalization condition introduced in ^7j) but a modification has been necessary to achieve 
good statistical precision. Section |3] contains the numerical determination of the scale depen- 
dence of the current in the SF scheme which is independent of the lattice formulation. We 
also relate the current renormalized at some proper low scale to the RGI current. Section 
gives our results for the lattice formulation dependent values of the Z-factor at this low 
scale. In Section |S] we then discuss Eq. p.ljl and explain in detail how our results are to be 
used. As an example we obtain -Fp^ at from published numbers of the bare matrix element. 
We finish with a brief discussion of the results in Section |HI Details of the numerical and 
perturbative calculations are described in appendices. 
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2 Intermediate renormalization scheme 



In this section we introduce our intermediate renormalization scheme. For reasons to be 
explained below, it differs from the one originally introduced in Ref. JJj. The perturbative 
calculations of Ref. are extended to the new scheme in Appendix B| 



We choose a mass-independent renormalization scheme, leading to simple renormaliza- 
tion group equations. The scheme is defined using the Schrodinger functional (SF) |18U19j . 
i.e. the QCD partition function Z = f TxL3 DL4, ip, ip] e - s i A ^^] Q n a T x L 3 cylinder in 
Euclidean space, where periodic boundary conditions in the spatial directions of length L 
and Dirichlet boundary conditions at times xq = 0, T are imposed on the gluon and quark 
fields. 1 Their exact form can be found in Ref. [25- Moreover, we set T = L throughout, 
and the renormalization scale \i is identified with the inverse box extension, 1/L. Such a 
finite- volume renormalization scheme is chosen, since it allows to study the scale dependence 
in the continuum limit for a large range of /i [22 ,23,24 , 13 . We can then relate the quantities 
renormalized at some low scale \i to the RGI quantities. Reviews of the strategy are found 
in |25U26llT4"| . and for a more detailed description the reader should consult Ref. ^3] which 
we will follow quite closely. 

As detailed in ^T%\, we consider the SF with vanishing boundary gauge fields and 9 = 0.5. 
These settings are identical to those used for the quark mass renormalization in [Oj and were 
motivated by meeting the criteria of good statistical precision of the Monte Carlo results, well- 
behaved perturbative expansions of the renormalization group functions and minimization 
of lattice artifacts [21] ■ Static quarks are included as discussed in Ref. |17j . and we use the 
notation of that reference. Throughout most of this paper, we formally stay in the framework 
of continuum QCD; some notation and basic formulae of the lattice regularized theory, in 



which the following expressions receive a precise meaning, are collected in Appendix A 

In contrast to the relativistic current, there is no axial Ward identity which protects the 
renormalized static-light axial vector current, 

(A$*) (x) = n tat ?i(^)7075^h(x) , (2.1) 

from a scale dependence. Its scale evolution is governed by the renormalization group equa- 
tion 

M^j=7(S)*, (2-2) 

where 

$ = $ stat = (f|(^ at ) |i) (2.3) 

is an arbitrary matrix element of the renormalized static current. The renormalization group 
function 7, the anomalous dimension, has a perturbative expansion 

7(5) ~t { 70 + lif + 725 4 + 0(5 6 ) } (2.4) 
with a universal leading coefficient |28[ I29| . 

70 = "i' (2 - 5) 

1 The spatial boundary conditions of the quark fields are only periodic up to a global phase 8 $13; an 
additional 'kmematical' parameter. 
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and 71, 72, . . . depending on the chosen renormalization scheme. 

Non-perturbatively, one computes the change of $ under finite changes of the renormal- 
ization scale. For a scale factor of two, the induced change defines the step scaling function, 

at\u) = *W2)/*(M) = Zr\2L)/Zr\L) , (2.6) 

whose argument u = g 2 (L) is taken to be exactly the coupling defined in |24| . and as always 
in the SF we have fi = 1/L. 



2.1 The old scheme 

In Ref. , a normalization condition was formulated in terms of suitable correlation func- 
tions defined in the SF. It reads 

Zf at (L) X(L) = A (0) (L) at vanishing quark mass , (2.7) 

with 

*M = S ^ L < 2 - 8 > 

and X(°)(L) the tree-level value of X{L). Here, /X tat is a correlation function between 
a static-light pseudoscalar boundary source and Ag tat in the bulk, and ff 1 ^ denotes the 
correlator between two such boundary sources at xq = and xq = T: 

ft\xo) = -^|d 3 yd 3 z<Ar t (x)C h (y)75Ci( Z )>, (2.9) 
/f at = ~ J d 3 ud 3 vd 3 yd 3 z <Ci / (u) 75 Ch / (v)Ch(y)75Ci(z)> . (2.10) 

(For the proper definition of the 'boundary quark and antiquark fields' £, £ we refer to 
Refs. |17|l21j .) The two correlators are schematically depicted in Figure and their explicit 
form on the lattice is given in Eqs. {QUI and (TP3|) of AppendixEj] In the ratio (JTHJ) both 




Figure 1: Sketch of the correlation functions f^ at (xo) (left) and /^ tat (right). The single 
and double lines represent light (i.e. relativistic) and static quark propagators, respectively. 

the multiplicative renormalization of the boundary quark fields and the mass counterterm 
of the static field cancel. 

We now have to point out a drawback of this scheme that only becomes evident, 
when it is implemented numerically. Namely, the lattice step scaling function X^ at (u, a/L) 
(cf. Eq. ()A.20|) for its definition), computed by means of Monte Carlo simulations, has large 
statistical errors at u ~ 1.5 and larger. In particular, these errors grow with L/a. This can 
be inferred from the results tabulated in Appendix IA.2I and is illustrated for three repre- 
sentative coupling values in Figure^ For L/a = 12,16 (which amounts to also calculate 
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Figure 2: Comparison of the numerical precision of the lattice step scaling function com- 
puted in the scheme of Ref. \171j with the new one introduced in Subsection \2.3l which will 
be used in the rest of this paper. (The symbols are slightly displaced for better readability) 
The statistics of both computations is the same. 



^stat £ Qr 2L/a = 24,32), already around u ~ 1.5 a precise determination of X^ at (ii, a/L) 
with a reasonable computational effort becomes impossible. The reason for this lies in the 
boundary correlator ff ta,t being part of the renormalization condition Eq. I)2.7j) : it contains 
the static quark propagating over a distance T = L. Thus J| tat falls roughly like e~ ei9 o T / a 
with ei = yr^p- x 19.95 j3U] the leading coefficient of the self-energy of a static quark. On 
the other hand, the statistical errors fall much more slowly, leaving us with an exponential 
degradation of the signal-to-noise ratio. 

To circumvent this problem, we now introduce a slightly modified renormalization 
scheme. (Therefore, for the rest of the paper, the scheme of Ref. ^7j - - if mentioned at 
all — will only be referred to by labelling the corresponding quantities with an additional 
'old', e.g. Zf at — * ZaoW-) The general idea is to replace /f tat containing a static and a light 
quark by two boundary-to-boundary correlation functions. One of them contains a light 
quark- antiquark pair, the other a static quark-antiquark pair. Both can be computed with 
small statistical errors, the latter because the variance reduction method of Ref. |31| can be 
applied. Since the static-static boundary correlator has not been studied before, we discuss 
it in some detail. 



2.2 The static-static boundary correlator /J 1 

We define 

/l h (*3) = -^2 / dx ± d *2 d 3 y d 3 z (Ch'(x)75Ch'(0) Ch(y)75Ch(z)> , (2.11) 

represented graphically in the left part of Figure G3 After integrating out the static quark 
fields, fi^ becomes the (trace of the) product of a timelike Wilson line and the complex 
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Figure 3: The correlation functions fi(xs) (left) and f\ (right). The notation is the same 
as in Figure 01 

conjugate one from boundary to boundary (see Eq. HB,1|) ). They are separated by x in space. 
This quantity is integrated over x\,X2 but retains a dependence on X3. In the following we 
take d = \xs\ as its argument, where the periodicity of the system in the space directions 
restricts it to < d < L/2. 

Upon renormalization the correlation function f± h becomes 

(/ 1 hh ) R (d) = e~ 2SmxL (Z h ) 4 f^(d) , 

where is the wave function renormalization constant of a static boundary quark field and 
5m the linearly divergent static mass counterterm. Therefore, to study the properties of 
further, we form the finite ratio 

,,,,, (/ri»w /rw t - 2 , n ,,,,, 

hW L,„) = m - m = 1Fm at ,<L)=„. (2.12) 

Considered on the lattice, it has a continuum limit. As outlined in Appendix IB.ll the 
one-loop coefficient h^(d/L) of the perturbative expansion 

h(d/L, u) = 1 + u (d/L) + u 2 hW (d/L) + ... (2.13) 

is given by 

-§(!-£)'. (2,4) 

Remarkably, this form holds exactly on the lattice without any a/L-dependence. Some 
insight why this is so is presented in the appendix as well. At two-loop accuracy, we do not 
expect exact a-independence any more, but still one may hope that the favourable kinematics 
keep lattice artifacts small. 

The one-loop expression is compared to results from our non-perturbative computation 
of for two representative values of the coupling in Figure |IJ The figure contains non- 
perturbative results for L/a £ {12,16,20,24,32} but at the level of our statistical errors, 
which are about 1% and smaller, no lattice artifacts of the ratio h can be seen. 

For low d, the non-perturbative data for h are well described by cx (L/d) + c' , where the 
constant c grows with u = g 2 . Hence the correlation function contains a non-integrable short- 
distance singularity, which is the reason why we will not integrate over d in the following. 
It is easy to see that this singularity is absent up to and including the order u 2 , but in 
higher-order terms in perturbation theory such a singularity may appear. 
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Figure 4: Non-perturbative d-dependence of h(d/L,u) = /| lh (d)//| lh (L/2) at two (low and 
high) values of the coupling, compared to one-loop perturbation theory. As data points 
corresponding to various lattice resolutions are included, the figure also reflects the weak 
cutoff dependence of this ratio. 



2.3 The new renormalization scheme 

Choosing d at its maximum value to further keep discretization errors at a minimum, we 
specify our (non-perturbative) renormalization scheme by 

Zf at (L) E(L) = H^(L) at vanishing quark mass . 

with 



(2.15) 



f^(L/2) 



[hf^(L/2)] 1/A ' 

Here, /i is the correlator between two light-quark pseudoscalar boundary sources, 
H = / d 3 ud 3 vd 3 yd 3 z (Ci / (u)75C2 / (v)C 2 (y)75Ci(z)} , 



(2.16) 



(2.17) 



depicted in the right part of Figure The form of fx and /J 1 * 1 on the lattice is given in 
Eqs. (|A.12|) and (|A.14|) of Appendix IA.11 As before, the combination of these correlators 
in the denominator of (|2.16j) is such that the boundary field renormalizations and the mass 
counterterm drop out and no other scale but L appears. 

For 9 = 0.5 the perturbative calculation summarized in |Appendix B now yields 

-^2 { 0.10(2) - 0.0477(13)JV f } , (2.18) 



7f 



(4vr) 5 



which differs only little from the one in the old scheme |17j . 

Note that O(a) improvement |32U21| can be applied and is an important ingredient in 



practice to reduce the cutoff effects in the numerical simulations (see Appendix A ) . Return- 
ing to Figure El one observes that the statistical errors of the lattice results are indeed much 
smaller in the new scheme. 
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3 Non-perturbative running and renormalization group invariant 



In this section we present our quenched results on the evolution of over more than 

two orders of magnitude in [i. To this end we consider the evolution of Zj£ under repeated 
changes of the scale (i.e. the box size L) by a factor of two at fixed bare parameters. Starting 
at some initial low-energy value (i.e. some large L = L max ), one thereby climbs up the energy 
scale by repeated application of the inverse of the step scaling function until the perturbative 
domain at high energies (i.e. small 2 _fc L max ) is reached, where finally the associated (scale 
and scheme independent) renormalization group invariant may be extracted. As in the 
previous section we keep the discussion in the continuum theory here; the underlying lattice 



calculations are described in Appendix A 



3.1 Step scaling function 



0.99 



0.96 

3 

"> 0.93 

3 



0.9 



0.87 



1 

- i — 


i i i 1 i 




1 1 1 


1 1 1 1 
u=1.0989 _ 






— a---- 




CD 




— s $ 






u= 1.5553 

03 _ 


^ ' EE H — _ ?R ii-O AAQA 










3 










.. _ ii = 2.770 


" 1 


i i i 1 i 


i i 


, 1 , 


u=3.480 

i - 

i i --..I i 





0.01 




0.02 


0.03 



(a/L)« 

Figure 5: Lattice step scaling function £^ at and its continuum limit extrapolations for 
some selected values of u. 

The evolution of Z^ at from size L to 2L is given by its step scaling function, <T^ at («), 
which has already been introduced in Eq. ([2.6)1 . but where it is understood that is 
defined in the new renormalization scheme according to Eq. ()2,15|) . 

As detailed in Appendix lA.21 the sets of lattice parameters (L/a, (3, k), which in practice 
are required to non-perturbatively compute cr^ at (u), can be taken over from the quark mass 
renormalization ^3]- The available coupling values u allow to trace the scale dependence of 
r^stat U p to L = 2L max , where the scale L max is implicitly defined through 



5 2 (£max) =3.48. 



(3.1) 



The sequence 



u k = g 2 (2 k L rs 



k = 0, . . . 



(3.2) 



7 



is known from Ref. ^3], and thus the corresponding sequence 

^tat^-fc+l^) $( 2 fc-l/L max ) 



i 



1, 



is simply given by 



v = l, v k+1 



a 



stat 



(Uk) 



(3.3) 
(3.4) 



once the function cr^ at (it) is available in the corresponding range of u. 



u 


cr^ at (it) 


1.0989 


0.9796(20) 


1.3293 


0.9746(25) 


1.4300 


0.9719(26) 


1.5553 


0.9668(27) 


1.6950 


0.9727(28) 


1.8811 


0.9632(28) 


2.1000 


0.9589(35) 


2.4484 


0.9432(27) 


2.7700 


0.9423(41) 


3.4800 


0.9154(41) 



Table 1: Results for the continuum step scaling function o"^ at (it) 
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Figure 6: Continuum step scaling function crj at (u) and its polynomial fit. 

The calculation of the lattice step scaling function and its subsequent continuum ex- 
trapolation yields the pairs u and <7^ at (it) listed in Table ^ An impression of the quality 
of the continuum extrapolation is gained from Figure El but for a more detailed account 
of the lattice simulations and data analysis we refer to Appendix IA.2I An interpolating fit 
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of a s ^ at (u) is shown in Figure H3 The leading coefficients (so and sx, see Appendix IC.1|) of 
the interpolating polynomial are fixed to the perturbative predictions, Eqs. (|B.14jl . This fit 
is then inserted into the aforementioned recursion, and propagating all errors through the 
recursion we obtain 

Zl tat (2L max )/Zl tat (L) = 0.7551(47) at L = 2~ 6 L max , (3.5) 

with the value of the coupling at this box size being g 2 (2~ 6 L max ) = 1.053(12) ^3]- Let us 
emphasize once more that L = 2L max and L = 2 _6 L max represent low- and high-energy 
scales, respectively, which in this way have been connected non-perturbatively. (Our data 
actually allow to go up to L = 2~ 8 L max .) 



3.2 RGI matrix elements of the static axial current 

We now proceed to relate the renormalized matrix element 

$Gu) = Zl tat (L)1>ba r c(9o) , H = l/L, 
at L = 2L max to the renormalization group invariant one defined by 2 



$rgi = x [2b g 2 (n) 



-70/2&0 



exp 



dg 



7(g) _ _To_ 
I3{g) b g 



with the universal leading-order coefficients bo = ll/(47r) 2 and 70 = — l/(4-7r 2 ) of 
7-functions, respectively. Casting this equation in the form 



(3.6) 

(3.7) 
the (5- and 



RGI 



$((2L max )-!) 



Zj^ at (2L max ) 

x [26 05 2 (/x)]- 70/2b ° 



exp 



9(m) 



dg 



7(ff) _ 7o_ 
/5(g) b g 



with fj, = 2 6 /L 

max! we see that the first factor is known from Eq. (|3,5j) . while in the second one 
only couplings g 2 < 1.05 contribute and it can safely be evaluated by perturbation theory. 
Still, for the perturbative error to be negligible, 7 has to be known to two-loop accuracy 
and (5 to three-loop. Upon inserting g 2 (2 6 /L max ) = 1.053 and numerical integration of the 
second factor we find 



$(/x)/$rgi = L 



at fi = (2L E 



(3.8) 



in the SF scheme. Entirely consistent numbers, with slightly larger errors, are obtained for 
<&(//) /<3?rgi if one switches to perturbation theory at fi = 2 7 /L max or fi = 2 8 /L max instead. 

In Figure we compare the numerically computed running with the corresponding 
curves in perturbation theory. While good agreement with the perturbative approximation 
is seen at high scales, a growing difference of up to 5% becomes visible when [i is lowered to 
fj, ps 2.5A. 

Below it will be more convenient to specify the scale fi in Eq. (|3.8|) in terms of tq |33j 
instead of L m3iX . Taking also the small error contribution from the uncertainty of -L max in 



In a loose notation, we take sometimes L and sometimes fi = l/L as the argument of g. 
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Figure 7: Numerically computed values of the running matrix element of the static axial 
current in the SF scheme compared to perturbation theory. The dotted and solid lines are 
obtained from Eq. i)^.7|) using the 1/2- and 2/3-loop expressions for the 7- and (3-functions, 
respectively, as well as AL max = 0.211 from Ref. \13tj . 

units of r , L ms , x /r = 0.718(16) [31], into account, the final result for the regularization 
independent part 3>(/i)/3>RGi of the total renormalization factor is 



Note that this result refers to the continuum limit so that the error on 3>(//)/$rgi oi about 
0.9% should only be added in quadrature to the proper matrix element under study after its 
continuum extrapolation. 

4 Zj^ at at low scale and total renormalization factor 

We still need to relate (^4f£ at )o(//), renormalized at some appropriate scale fi, to the bare 
lattice operator. This amounts to computing Z^ at at the low-energy matching scale L = 
2-^max = 1.436 ro, which is briefly explained in Appendix IC.2I Since in this step the bare 
operator is involved, Z^ at does depend — in contrast to the result of the previous section 
- on the choice of action. We have considered three different cases. The first two are the 
non-perturbatively O(a) improved action of Ref. with c^ at = — x (= one-loop) and 
separately with c^ at = 0. Their combination will in the future allow to study the influence 
of c^ at on the continuum extrapolations of renormalized matrix elements. The third choice 
is the unimproved Wilson action which is of interest, because so far the best computations 
of the bare matrix element did not use improvement |101lllj . 

The numerical results for Zj£ are shown in Figure El For later use they are represented 
by interpolating polynomials, 



*(M)/$RGi = 1-088(10) at fi = (1.436 r ) 



1 



(3.9) 




(4.1) 
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with coefficients Z{ as listed in Tabled The statistical uncertainty to be taken into account 
when using this formula is about 0.4%. 






_stat 
C A 


applicability 


i 


Zi 


fi 


non-perturbative 


-5oV(4vr) 


6.0 < < 6.5 





0.6944 


0.6382 


m 






1 


-0.0946 


-0.0869 








2 


0.1239 


0.1139 


non-perturbative 





6.0 < < 6.5 





0.6750 


0.6204 








1 


-0.0838 


-0.0771 








2 


0.1200 


0.1103 


0, 





5.7 < < 6.5 





0.6090 


0.5598 








1 


-0.1186 


-0.1090 








2 


0.3438 


0.3160 








3 


-0.2950 


-0.2711 



Table 2: Coefficients of the interpolating polynomials, Eq. i)4.1j) and Eg. i|4.<jj) . Uncertainties 
are discussed in the text. 



The total renormalization factor 



The total renormalization factor to directly translate any bare matrix element $bare(5o) of 
^stat ^ n ^ j^qj m atrix element, $rgi> can be written as 



-Zrgi(S'o) 



RGI 



^=(1.436 ro)" 



x Zf a \g ,L/a) 



L=1.436 r 



(4.2) 
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We combine Eq. (|4.1|) with Eq. Q3.9I) and represent the total Z-factor by further interpolating 
polynomials, 

^RGi(5o) = ^/ i (/?-6) i , (4.3) 

i>0 

whose coefficients are also found in Table |2 These parametrizations of Zrgi are to be 
used with an uncertainty of about 0.4% 3 at each /3-value and an additional error of 0.9% 
(from <&hgi/&(h)), which remains to be added in quadrature after performing a continuum 
extrapolation. 

5 Matrix elements at finite values of the quark mass 

In order to use results from the static theory, one still has to relate its renormalization group 
invariant matrix elements to those in QCD at finite values of the quark mass, to. This step 
may also be seen as a translation to another scheme, defined by the condition that matrix 
elements in the static effective theory renormalized in this scheme and at scale fi = to are 
the same as those in QCD up to l/m-corrections. This scheme is therefore denoted as the 
'matching scheme' )14| . Below, we will specify precisely which quark mass to is to be taken. 

5.1 Conversion to the matching scheme 

Let us write the relations for the special case of the matrix element of the axial current 
between the vacuum and the heavy-light pseudoscalar, 

$RGI = ^RGl(PS|< at |0>. (5.1) 

We then have 

FpsV^vs = Cp S (to) x $ RG i + 0(1/to), (5.2) 

where to is the MS quark mass at renormalization scale to. 4 The function Cps(^) is given 
by 



5ps(/,) = [26o^(^)] 7o/2bo exp 




7(g) _ 7o_ 
0{g) b g 



(5.3) 



with g(/i) the MS running coupling and 7 the anomalous dimension in the matching scheme. 
The latter is known to two loops [2211123111302]] with 70 being the same as before and 

71 EE 7r tch =7 MS_i5_ ) ( g_ 4) 



,ms 1 ' 



^ - -57^(- +28 « 2) - 5iV 'J' (5 ' 5 > 

For illustration, Cps(/^) is plotted in the upper part of Figure where for the numerical 
evaluation the /^-function is always taken at four-loop precision jl^, while to estimate the 

Only in the C3,SG Csw = the error to be associated with the formulae for Z^ at and Zb,gi grows to 0.5% 

at P « 6.3 and 0.8% at /3 a 6.5. 

4 Note that in |36| a similar equation with the pole mass instead of the MS mass is written. At the two-loop 
order, which will be used below, this does formally not make any difference. However, the pole mass does 
not have a well-behaved perturbative expansion |37|. and we therefore prefer a short-distance mass such as 
the MS mass. 
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perturbative uncertainty we show the result for the one-loop and the two-loop approximation 
of 7. 



Cp S (m) 



1.4 
1.3 
1.2 
1.1 h 




_l I I I I I I I l_ 



_L 



0.2 0.4 



0.6 



0.8 m- 



C PS (A^/M) 



[GeV" 1 ] 



1 1 1 1 1 1 1 1 r 



1- loop 7 

2— loop 7 




Figure 9: Conversion factor to the matching scheme, which translates the RGI matrix 
element to the one at finite mass. 

Eq. 1)5.3)) can be rewritten in a form displaying explicitly that also this step is not 
restricted to perturbation theory. In terms of the renormalization group invariant quark 
mass, M, we have 



FpsV^ = C PS (M/A m ) x $ RGI + 0(1/M) , 



(5.6) 



where now only renormalization group invariants enter. To evaluate Cps(M/Aj^g) in pertur- 
bative approximation, one changes the argument of the function Cps, Eq. (|5.3jl . by inserting 



M/m^) = [2b g 2 (ji) 



-do/2&o 



cxp 



dg 



dp 
bog 



(5.7) 



together with the condition m(m^) = where r(g) denotes as usual the renormaliza- 

tion group function of the renormalized (running) quark mass with universal leading-order 
coefficient do = 8/(47r) 2 . A numerical evaluation (with the four-loop r-function |42U43j in 
Eq. 1)5.71) ) is shown in the lower part of Figure E3 Eq. (|5,6|) is not only the cleanest form 
from a theoretical point of view but it is also practical, because the relation between bare 
quark masses in the 0(a) improved lattice theory and the renormalization group invariant 
mass M is known non-perturbatively in the quenched approximation (441145) . 

For later convenience we represent Cps in terms of the variable x = 1/ln (M/A^jg) in 
a functional form motivated by Eq. (|5.3|) . 



C PS = x lo/2bo { 1 - 0.065 x + 0.048a; 2 } , x = 1/ln (M/A m ) < 0.52 , 



(5.8) 
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with bo = ll/(47r) 2 and 70 = — l/(4-/r 2 ). It describes the result for the two-loop approxima- 
tion of the 7-function within less than 0.01%. Of course in this step a perturbative error is 
involved, which is difficult to estimate. Assuming a geometric growth of the coefficients of 
the 7-function, we find that the 72-term would cause a change by around 1% at m = m b yg 
and by 2.5% at m = 1.2 GeV. Thus one may attribute a 2-4% error due to the perturbative 
approximation, which could be much reduced by a computation of the three-loop anomalous 
dimension. 

In principle, Cps may be computed also non-perturbatively following the strategy out- 
lined in Ref. |14j . It is then defined only up to 1/M-terms, consistent with Eq. (|5.6j) . 

5.2 Application: first non-perturbative renormalization of -Fg* at 

We now take bare matrix elements of ^4q for unimproved Wilson fermions from the litera- 
ture to obtain an estimate for F^ a in the static approximation. This exercise serves mainly 
to illustrate how to use our results. 

1. The matrix elements are needed at a fixed value of the light quark mass. To avoid 
issues in the extrapolation to very light quarks, we here consider only F^ s . To fix the strange 
quark mass, we use that the sum of the light quark masses is to a good approximation 
proportional to the squared (light-light) pseudoscalar masses, mp s (l, 1) [31], and interpolate 
the data for the decay constant of |10[lll| as a function of mp s (l, l)rQ to mp s (s, s)rg = 
(2m\ - m\)rl = 2m\rH(\ + m x /m s ) = 3.0233. (To arrive at the latter, we employed 
m K r o = 1-5736 |44| and m s /m\ = 24.4 from chiral perturbation theory |46|.^ The resulting 
dimensionless numbers rp^^bare are listed in Table El 



P = 6/5o 2 


5.7 


5.9 


6.0 


6.1 


6.3 


ro/a 

r0 /2 $bare 


2.93(1) 

4.75(25) 

2.99(16) 


4.48(2) 

4.09(21) 

2.35(12) 


5.37(2) 

3.94(13) 

2.21(7) 


6.32(3) 

3.79(36) 

2.09(20) 


8.49(4) 

4.00(29) 

2.20(16) 



Table 3: Matrix elements of ylg tat m units of ro for the light quark mass equal to the strange 
quark mass. Bare matrix elements come from Ref. \lUf . with the exception of [3 = 6.0 which 
is taken from Ref. The scale ro/a 1331 is used as determined in Ref. \34\ . and its 

uncertainty is already included here. 

2. We renormalize by multiplying with Eq. (|4.3|) . using the /, from Table |U(c sw = ca = 
gstat _ o) 5 take into account a 0.4-0.5% error from the non-universal part of the Z-factor 
at each value of go and find the last line in Table |21 Assuming the leading linear behaviour 
in the lattice spacing to dominate for a/ro < 1/4, we extrapolate to the continuum limit as 
shown in Figure 1101 Adding to the extrapolation error in quadrature also the 0.9% error 
contribution of ZrgIj which is independent of go, yields 

r-o /2 $RGi = 1-93(34) at a = 0. (5.9) 

3. Finally, inserting M h r = 17.6(5) gHHl] and A^ro = 0.602(48) 13,, one gets via 
the formula in Eq. (|5.8j) 

C PS (M b /A M s) = 1.23(3), (5.10) 
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Figure 10: Continuum extrapolation of the non-perturbatively renormalized matrix element 
of Af at based on the unimproved Wilson data for F^ are from Refs. \1 (AM 1] . 

where a 2% error for the perturbative approximation is assumed. With the experimental 
spin-averaged B-meson mass me = w-b s = 5.4 GeV, we then obtain from Eq. (|5.6|) : 

F£ at r = 0.64(11), (5.11) 
F B tat = 253(45) MeV for r = 0.5fm. (5.12) 

The result contains all errors apart from the uncertainty owing to the quenched approxima- 
tion. Evidently, the continuum extrapolation may be done much better, once O(a) improved 
results with sufficient precision and small lattice spacings are available. 

6 Discussion 

We have performed the scale dependent renormalization of Aq 1 ^ by constructing a non- 
perturbative renormalization group in the Schrodinger functional scheme, and agreement 
with perturbation theory at large scales was demonstrated. The renormalization factors 
needed to extract the associated RG invariant are computed with good numerical accuracy, 
which is a crucial prerequisite for a controlled determination of F-q in the static limit. In 
Ref. ^3] it was shown that the renormalization factors obtained in this way differ appreciably 
from earlier estimates ^U] based on tadpole-improved perturbation theory jlHj. Hence their 
non-perturbative computation is important. 

We have not emphasized this so far, but our computation provides the scale dependence 
of all static-light bilinears 

Or(x) =^ l {x)TMx) ■ (6.1) 

They are renormalized by 

(0 7 J R (z) = Zg*0 lh , 

(0 7o ) R (z) = ZZl tat O 70 , (0 7fc75 ) R (x) = ZZl tat 7 , 75 , (6.2) 

with a scale independent renormalization Z. This pattern is due to the heavy quark spin 
symmetry which is exact on the lattice, and due to the chiral symmetry of the continuum 



t | i i i i | i i i i | i i i i | i i i r 
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theory. The latter means that the relative renormalization Z may be fixed by imposing a 
suitable chiral Ward identity 49 and is thus scale independent. 

Returning to the case of most interest, i ? g tat , our continuum extrapolation in Section T5. 21 
that uses unimproved data for the bare matrix elements from the literature and also quite 
large lattice spacings leaves much room for improvement of the present result, i*g tat = 
253(45) MeV. Apart from the obvious step of obtaining 0(a) improved bare matrix elements 
at small lattice spacings and extrapolating to the continuum, it will be necessary to estimate 
the 0(1/M) correction. There are two possible roads towards this goal. 

An elegant and clean way is to compute the 1/M-corrections directly as perturbations 
to the static effective theory. Again, the main problem here is renormalization. Indeed, 
this is a severe one, since mixing between operators of different dimensions has to be taken 
into account. This will require much more theoretical and numerical effort; but a possible 
strategy exists |i7| H4*]. 

In the mean time, one may also compare the prediction from the static approximation 
to what one obtains at M ~ M c and also to the results obtained directly at M = Mb, 
most notably the ones of Refs. [2J|H]. As emphasized in Refs. [57 H l5 H l52| this should be done 
in the continuum limit, since 0(a) errors get enhanced when the quark masses increase. 
At the charm quark mass these are sizeable but can be extrapolated away, at least in the 
quenched approximation [53] . The comparison between finite-mass decay constants and 
Fg tat is most conveniently done by comparing Fpsy / mps/Cps(M/Aj^g). Unfortunately, at 
present the error of the static result is still too large to draw any strong conclusions about 
1/M-corrections in F B : 

rl /2 $rgi = 1.93(34) , static : Eq. Q , (6.3) 
3/2 F BBy /m^ 



C PS (M b /A- 



1.46(23), using F Bs = 192(30) MeV 0E] , (6.4) 



r o /2 7^%T7?^T = L29 ( 5 )> usin § ^D S = 252( 9) MeV M- (6-5) 

Still, the difference of Eq. ()6.3|) and Eq. ()6.5|) shows that there are significant 1/M-corrections 
in the charm mass region. 

As a more technical remark we point out that the function h(d/L,u), Eq. (|2.12|) . shows 
very small a-effects in the quenched approximation and may be worth studying to verify 
improvement with dynamical fermions. 



Acknowledgements 

We thank DESY for allocating computer time on the APE-Quadrics computers at DESY 
Zeuthen to this project. This work is also supported in part by the EU IHP Network on 
Hadron Phenomenology from Lattice QCD under grant HPRN-CT-2000-00145 and by the 
DFG Sonderforschungsbereich SFB/TR 9. 



Appendix A Computation of the lattice step scaling function 

This appendix describes some details of the numerical simulations on the lattice and the 
subsequent calculations that we have performed in order to determine the step scaling func- 
tion for Z| tat . At the beginning we also recall a few basic definitions and formulae, which 
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are specific to the inclusion of static quarks and the correlation functions that are considered 
in the framework of the Schrodinger functional (SF). As the impact of the static quarks on 
0(a) improvement of the static-light sector has extensively been discussed in Ref. the 
reader might consult this reference for further details and any unexplained notation. 



A.l Definitions 

Lattice action 

The total lattice action is given by the sum 

SlUrf^faM =S G [U]+S F [U,^ h ^}+S h [U,^ h ,xlj h }, (A.l) 

where 5g and Sf are the standard pure gauge and 0(a) improved Wilson actions for rela- 
tivistic (light) quarks, see Eqs. (A. 22) and (A. 23) - (A. 26) of Ref. respectively, and Sh 
denotes the lattice action for the heavy quark: 

S h [U3 h ,^h] = a 4 ]T^ h (x)VuVh(x) • (A.2) 

X 

The fields Vh an d V'h °f the static effective theory are constrained in such a way (namely 
P+tph = V'h and ip^P+ = "0h) that one is left with just two degrees of freedom per space-time 
point 4 a and only the (time component of the) backward lattice derivative, V^, enters in 
the action (jA.2|) , Hence, static quarks propagate only forward in time, which also reflects in 
the form of the associated quark propagator, 

S h (x,y) = U(x-a6,Oy 1 U(x-2a6,Oy 1 ■■■ C/(y,0) _1 

x 8(x Q - y )5(x - y)P+ , P± = |(1 ± 7o ) , (A.3) 

being just a straight timelike Wilson line. 

To impose SF boundary conditions, Eqs. (|A.1|) and ()A.2|) are supplemented by 

ipi(x) = if xo < or xo > L , 

iph(x) = if x < or x > L (A. 4) 



and 



P-ik{x) =P+Mx) =0, (A.5) 



XQ = L 



10=0 

while in the pure gauge part the spatial plaquettes at xo = and xq = L receive a non-trivial 
(and coupling dependent) weight, see Eq. (A. 29) of Ref. |T7j. In the discretization of the 
SF as described in we choose zero boundary gauge fields throughout, C = C = 0, 

which translates into the boundary conditions U(x, k)\ xo= o = U(x, k)\ Xo= L = 1 for the lattice 
gauge field. Similarly, the light and static quark fields at xq = 0, L are fixed to appropriate 
(space dependent) boundary functions; the corresponding boundary conditions are collected 
in Eqs. (3.2), (3.3) and (3.5) of Ref. JZj and not repeated here. 
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SF correlation functions 

Observables are then denned as usual through a path integral involving the total action S. In 
this work we focus on SF correlation functions that are constructed from the 0(a) improved 
static-light axial current 

(Af tat ) (a;) = ^ tat (x) + acT t 5A s t! "\x), (A.6) 
Af at (x) = Mxho^Mx), (A.7) 
5Af*\x) = ^( x ) 7i75 I(^. + V*)^h(x). (A.8) 

Unless it is indicated differently, the improvement coefficient c^ at is set to its one-loop 
perturbative value, 

«r = -^5o 2 , (A.9) 

computed in Refs. On the lattice, in terms of the boundary quark fields the 

correlation functions of these fields, as well as the various types of pseudoscalar correlators 
from one boundary to the other that are needed in addition, read explicitly: 



ft\x ) = -a 6 ^i(^rWC h (y)75Ci( Z )> , (A.10) 
fff(xo) = -a 6 ^i<o< at (x)Ch(y)75Ci(z)> , (A.ll) 



2 



A = E ^<C/(u) 7 5C/(v)C,(y)750( Z )> , (A.12) 

u,v,y,z 

/f at = E ^<Ci / (u)75Ch / (v)Ch(y)75Ci( Z )> , (A.13) 

u,v,y,z 

Zi hh (*3) = E ^<Ch'(x)75Ch / (0)Ch(y)75Ch(z)> . (A.14) 

Moreover we introduce the two ratios 

Xj(g , Lj a) -j== , (A.15) 

•=i(0o, La) = — —j- , (A.16) 

[fi fi (L/2)] 

which are constructed such that the (unknown) wave function renormalization factors of the 
boundary quark fields as well as the (linearly divergent) mass counterterm 5m cancel out 
and only the static current remains subject to renormalization. 

Renormalization 

In Ref. the renormalization constant Z^**" = entering the 0(a) improved static 

axial current renormalized in the SF scheme, 

(^g at )o = Zl tat (l + ^ tat am q )(^f at ) , (A.17) 



18 



was defined in terms of the ratio Eq. (jA. 151) by imposing the renormalization condition (with 
mo, m q and m c as defined in [T7j ) 

Zf% d (g , L/a) Xi(go, L/a) = Xj(0, L/a) at m = m c . (A. 18) 

Thus, naturally runs with the scale /i = 1/L. In the present context it will be referred 

to as the 'old' scheme, whereas the so-called 'new' scheme based on Eq. (|A.16|) is specified 
by 

Zt* t (g ,L/a)E 1 (g ,L/a) = E l (0,L/a), m = m c , L = 1/ 1 m. (A.19) 

For 6 = 0.5, which is chosen in our simulations, the relevant values of the tree- level nor- 
malization constant Ei(0, L/a) (or E^\a/L) in the notation of Appendix B summarizing 



the perturbative calculations) are collected in Table IA.11 As an aside we remark that 
E^O, L/a) = X l (0,L/a) holds. 



L/a 


Hi (0, L/a) 


6 


-1.5964837518021 


8 


-1.5996643156321 


10 


-1.6011462370857 


12 


-1.6019540566018 


16 


-1.6027594410020 


20 


-1.6031330222299 


24 


-1.6033361949926 


32 


-1.6035384024722 



Table A.l: The ratio Ei for 9 = 0.5 at tree-level. 

The critical quark mass is always understood to be defined from the non-perturbatively 
O(a) improved PCAC mass in the light quark sector as in Ref. (i.e. for 9 = and T = L, 
evaluating the associated combination of correlation functions at xq = T/2). 

Lattice step scaling function 

The lattice step scaling function of the static axial current is defined through 

S f(^^) = W^7T at 9 2 (L)=u,m = m c . (A.20) 
Z% ™ (g , L/a) 

The additional condition mo = m c from above, referring to lattice size L/a, defines the 
critical hopping parameter value, k = k c . Moreover, enforcing g 2 (L) to take some prescribed 
value u fixes the bare coupling value g^ = 6/f3 to be used for given L/a. In this way S^ at 
becomes a function of the renormalized coupling u, up to cutoff effects, and approaches its 
continuum limit as a/L — * for fixed u. 

A. 2 Simulation details and results 

As emphasized before, our quenched lattice simulation and the data analysis are analogous to 
Ref. except that the boundary coefficient ct is set to its two-loop perturbative value |5T| : 

c 2-ioo P = 1 _ Q Q89 g 2 _ o_030 gl . (A.21) 
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The boundary 0(a) improvement terms involving quark fields have to be multiplied with a 
coefficient ct, which is known to one-loop from |^H], viz. 

c t i-i°°p = 1 _ .018#§. (A.22) 

Of course, owing to a priori unknown precision to which perturbation theory approximates 
these coefficients, linear lattice spacing errors are not suppressed completely, and we will 
come back to this issue later. As for the other contributing 0(a) improvement coefficients, 
we used the non-perturbative values for c sw and ca of [3S] for the relativistic fermions and 
the one-loop estimate (|A.9jl for c^ at in the static- light axial current. 

The renormalization constants Z A tat (g , L/a) and Z A tat (g , 2L/a) in Eq. ()A.20|) have 
been evaluated from the correlation functions in Eqs. (|A.10|I - (|A.14j) . which were computed 
in a numerical simulation with 9 = 0.5. (The latter parameter specifies the boundary con- 
ditions of the quark fields, see e.g. [20].) These simulations were performed on the APE-100 
parallel computers with 128 to 512 nodes, employing for the updating of the gauge fields 
the same hybrid-overrelaxation algorithm as in [241113] with two overrelaxation sweeps per 
heatbath sweep within a full iteration. This mix of updating was found to be close to opti- 
mal in [50]. Since the computation of SF correlators has already been detailed in Refs. |35| 
and (Appendix A. 2. 2 of) |13| . we just mention that we differ from them only by using the 
implementation [SU] of the SSOR-preconditioned BiCGStab inverter jSJ to solve the lattice 
Dirac equation. 

The computation of /^(d), where d = \x$\ and a < d < L/2, amounts to evaluate 
Eq. (|B.1|) . In order to improve the statistical precision of fi, the links building up the 
observable are evaluated by a 10-hit multi-hit procedure [31] . where each hit consists of a 
Cabibbo-Marinari heatbath update in three SU(2)-subgroups of SU(3). Translation invari- 
ance is fully exploited. 

In order to keep autocorrelations small, the measurements of the correlation functions 
were always separated by L/(2a) update iterations (and, respectively, five iterations in the 
case of f^ h ). Statistical errors stem from a standard jackknife analysis, where we have 
checked explicitly for the statistical independence of the data by averaging them into bins of 
a few consecutive measurements beforehand. The total number of measurements itself was 
such that the statistical error of £ A at was dominated by the uncertainty in Z A tat (go, 2L/a). In 
general, the uncertainties in the coupling and in the value of k c would have to be propagated 
into the error of £ A at as well. But as the former can be estimated to be much smaller than 
the statistical error of E A at and E A at is found to depend rather weakly on the bare (light) 
quark mass, we neglected both contributions in the final error estimate. 

In Table lA~2l - Table lA~4l we list our results on the step scaling functions of the static axial 
current 5 and — since they are available from our computations as well — of the pseudoscalar 
density defined as in [T3"j . The values of f3 and the critical hopping parameter k = k c to 
be simulated were taken over from Ref. |T3] without changes, which means to stay with 
c t and ct to one- loop accuracy in realizing the conditions g 2 (L) = u and tuq = m c . Note 
once more, however, that here, in contrast to ^3], for the corresponding renormalization 
constants themselves — particularly when comparing the results for Zp and Sp quoted in 
that previous work with those of the present one — the two-loop formula for ct , Eq. (|A.21|) , 
has been used. 

5 Here we do not tabulate the results on the static-static boundary correlator /} lh separately, but the 
numbers can be obtained from the authors upon request. 
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?{L) 


P 


K 


L/a 


Zp(g ,L/a) 


Zp(g , 


2L/a) 


Ep(u, a/L) 


1.0989 


9.5030 


0.131514 


6 


0.8190 


10) 


0.7846 


;q) 


0.9581 


16) 




9.7500 


0.131312 


8 


0.8107 




0.7786 


12) 


0.9604 


18) 




10.0577 


0.131079 


12 


0.8012 




0.7702 


11) 


0.9613 


16) 




10.3419 


0.130876 


16 


0.7956 


:§) 


0.7625 


'!5) 


0.9584 


;21) 


1.3293 


8.6129 


0.132380 


6 


0.7930 


;io) 


0.7493 


10) 


0.9449 


17) 




8.8500 


0.132140 


8 


0.7827 


;n) 


0.7402 


11) 


0.9456 


19) 




9.1859 


0.131814 


12 


0.7737 


7) 


0.7353 


12) 


0.9503 


;!7) 




9.4381 


0.131589 


16 


0.7666 


;n) 


0.7274 


18) 


0.9488 


[27) 


1.4300 


8.5598 


0.132453 


8 


0.7702 


;io) 


0.7273 


13) 


0.9443 


[21) 




8.9003 


0.132095 


12 


0.7610 


[6) 


0.7223 


14) 


0.9491 


[21) 




9.1415 


0.131855 


16 


0.7555 


7) 


0.7123 


'20) 


0.9428 


[28) 


1.5553 


7.9993 


0.133118 


6 


0.7666 


7) 


0.7165 


16) 


0.9346 


[23) 




8.2500 


0.132821 


8 


0.7590 


:§) 


0.7134 


13) 


0.9399 


[20) 




8.5985 


0.132427 


12 


0.7473 


10) 


0.7035 


13) 


0.9414 


[21) 




8.8323 


0.132169 


16 


0.7421 


;9) 


0.6976 


19) 


0.9401 


[29) 


1.6950 


7.9741 


0.133179 


8 


0.7442 




0.6939 


;!5) 


0.9325 


[24) 




8.3218 


0.132756 


12 


0.7341 


7) 


0.6862 


'!5) 


0.9348 


[22) 




8.5479 


0.132485 


16 


0.7277 


13) 


0.6805 


18) 


0.9352 


[30) 


1.8811 


7.4082 


0.133961 


6 


0.7348 




0.6764 


;e) 


0.9205 


14) 




7.6547 


0.133632 


8 


0.7258 


7) 


0.6691 


15) 


0.9219 


[23) 




7.9993 


0.133159 


12 


0.7173 




0.6632 


;s) 


0.9245 


12) 




8.2415 


0.132847 


16 


0.7117 


13) 


0.6604 


;20) 


0.9279 


;33) 


2.1000 


7.3632 


0.134088 


8 


0.7088 


13) 


0.6433 


16) 


0.9076 


;28) 




7.6985 


0.133599 


12 


0.6971 


;s) 


0.6385 


;24) 


0.9160 


;36) 




7.9560 


0.133229 


16 


0.6919 


12) 


0.6303 


'!7) 


0.9110 


;29) 


2.4484 


6.7807 


0.134994 


6 


0.6845 


10) 


0.6110 


12) 


0.8925 


;2i) 




7.0197 


0.134639 


8 


0.6784 


;s) 


0.6061 


19) 


0.8933 


;30) 




7.2025 


0.134380 


10 


0.6733 


;s) 


0.6021 


12) 


0.8943 


;2i) 




7.3551 


0.134141 


12 


0.6722 


;n) 


0.6012 


12) 


0.8944 


[24) 




7.6101 


0.133729 


16 


0.6661 


;s) 


0.5962 


10) 


0.8950 


;!7) 


2.7700 


6.5512 


0.135327 


6 


0.6619 


10) 


0.5758 


;20) 


0.8699 


;33) 




6.7860 


0.135056 


8 


0.6541 


13) 


0.5751 


'!7) 


0.8792 


;3i) 




6.9720 


0.134770 


10 


0.6505 




0.5717 


'!7) 


0.8788 


[28) 




7.1190 


0.134513 


12 


0.6482 




0.5705 


10) 


0.8802 


19) 




7.3686 


0.134114 


16 


0.6442 


16) 


0.5668 


16) 


0.8798 


;33) 


3.4800 


6.2204 


0.135470 


6 


0.6173 




0.5067 


11) 


0.8208 


;2i) 




6.4527 


0.135543 


8 


0.6133 




0.5101 


;21) 


0.8316 


;35) 




6.6350 


0.135340 


10 


0.6112 


;n) 


0.5078 


19) 


0.8307 


;35) 




6.7750 


0.135121 


12 


0.6076 


7) 


0.5061 


14) 


0.8329 


[24) 




7.0203 


0.134707 


16 


0.6063 


7) 


0.5097 


11) 


0.8406 


[21) 



Table A. 2: Results for the step scaling function Sp. 
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f(L) 





K 


L/a 


Z^(g ,L/a) 


Z^(g ,2L/a) 




,a/L) 


1.0989 


9.5030 


0.131514 


6 


0.8926(7) 


0.8698(8) 


0.9745 


12) 




9.7500 


0.131312 


8 


0.8860(7) 


0.8668(11) 


0.9782 


14) 




10.0577 


0.131079 


12 


0.8800(6) 


0.8630(11) 


0.9806 


14) 




10.3419 


0.130876 


16 


0.8786(7) 


0.8586(16) 


0.9773 


[20) 


1.3293 


8.6129 


0.132380 


6 


0.8733(8) 


0.8458(9) 


0.9686 


13) 




8.8500 


0.132140 


8 


0.8677(9) 


0.8418(13) 


0.9701 


18) 




9.1859 


0.131814 


12 


0.8635(7) 


0.8401(13) 


0.9729 


^7) 




9.4381 


0.131589 


16 


0.8593(9) 


0.8361(19) 


0.9731 


[25) 


1.4300 


8.5598 


0.132453 


8 


0.8593(7) 


0.8328(13) 


0.9692 


;!7) 




8.9003 


0.132095 


12 


0.8545(7) 


0.8314(14) 


0.9731 


18) 




9.1415 


0.131855 


16 


0.8516(6) 


0.8238(23) 


0.9674 


;28) 


1.5553 


7.9993 


0.133118 


6 


0.8572(6) 


0.8246(13) 


0.9619 


16) 




8.2500 


0.132821 


8 


0.8517(6) 


0.8248(13) 


0.9684 






8.5985 


0.132427 


12 


0.8459(8) 


0.8190(14) 


0.9683 


;20) 




8.8323 


0.132169 


16 


0.8425(9) 


0.8141(20) 


0.9662 


;26) 


1.6950 


7.9741 


0.133179 


8 


0.8414(9) 


0.8069(18) 


0.9590 


;24) 




8.3218 


0.132756 


12 


0.8359(8) 


0.8074(13) 


0.9659 


19) 




8.5479 


0.132485 


16 


0.8329(13) 


0.8081(19) 


0.9703 


;28) 


1.8811 


7.4082 


0.133961 


6 


0.8362(7) 


0.7939(7) 


0.9495 


11) 




7.6547 


0.133632 


8 


0.8290(6) 


0.7903(17) 


0.9533 


[21) 




7.9993 


0.133159 


12 


0.8247(7) 


0.7907(12) 


0.9588 


16) 




8.2415 


0.132847 


16 


0.8221(13) 


0.7898(23) 


0.9607 


[2,2) 


2.1000 


7.3632 


0.134088 


8 


0.8193(10) 


0.7732(20) 


0.9436 


;26) 




7.6985 


0.133599 


12 


0.8117(9) 


0.7756(25) 


0.9555 


;32) 




7.9560 


0.133229 


16 


0.8081(12) 


0.7704(21) 


0.9533 


;29) 


2.4484 


6.7807 


0.134994 


6 


0.8035(8) 


0.7420(12) 


0.9235 


18) 




7.0197 


0.134639 


8 


0.7945(7) 


0.7444(19) 


0.9370 


;26) 




7.2025 


0.134380 


10 


0.7936(9) 


0.7431(17) 


0.9363 


;24) 




7.3551 


0.134141 


12 


0.7930(10) 


0.7465(17) 


0.9414 


;24) 




7.6101 


0.133729 


16 


0.7907(8) 


0.7444(16) 


0.9415 


[22) 


2.7700 


6.5512 


0.135327 


6 


0.7886(9) 


0.7133(19) 


0.9045 


;26) 




6.7860 


0.135056 


8 


0.7791(12) 


0.7187(24) 


0.9225 


;35) 




6.9720 


0.134770 


10 


0.7786(9) 


0.7223(23) 


0.9276 


;3i) 




7.1190 


0.134513 


12 


0.7740(11) 


0.7220(17) 


0.9329 


;25) 




7.3686 


0.134114 


16 


0.7755(16) 


0.7281(30) 


0.9388 


13) 


3.4800 


6.2204 


0.135470 


6 


0.7587(10) 


0.6562(39) 


0.8649 


;53) 




6.4527 


0.135543 


8 


0.7496(11) 


0.6624(29) 


0.8837 


11) 




6.6350 


0.135340 


10 


0.7477(11) 


0.6658(31) 


0.8906 


13) 




6.7750 


0.135121 


12 


0.7451(11) 


0.6696(22) 


0.8988 


;32) 




7.0203 


0.134707 


16 


0.7470(10) 


0.6792(24) 


0.9092 


;34) 



Table A. 3: Results for the step scaling function S^ at (in the 'new' scheme). 
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t(L) 


P 


K 


L/a 


Z&(9o,L/a) 


Z%% d (g ,2L/a) 


sSd(«.«/i) 


1.0989 


9.5030 
9.7500 


0.131514 
0.131312 


6 

8 


0.8903(8) 
0.8846(7) 


0.8676(8) 
0.8635(11) 


0.9745(12) 
0.9762(14) 




10.0577 


0.131079 


12 


0.8791(5) 


0.8628(14) 


0.9814(17) 




10.3419 


0.130876 


16 


0.8773(9) 


0.8601(21) 


0.9804(26) 


1.3293 


8.6129 


0.132380 


6 


0.8721(9) 


0.8430(10) 


0.9666(16) 




8.8500 


0.132140 


8 


0.8664(10) 


0.8428(14) 


0.9728(19) 




9.1859 


0.131814 


12 


0.8616(6) 


0.8351(19) 


0.9693(23) 




9.4381 


0.131589 


16 


0.8580(10) 


0.8370(48) 


0.9755(57) 


1.4300 


8.5598 


0.132453 


8 


0.8577(8) 


0.8325(14) 


0.9706(19) 




8.9003 


0.132095 


12 


0.8535(6) 


0.8332(21) 


0.9761(26) 




9.1415 


0.131855 


16 


0.8488(6) 


0.8172(55) 


0.9628(65) 


1.5553 


7.9993 


0.133118 


6 


0.8552(7) 


0.8243(19) 


0.9639(23) 




8.2500 


0.132821 


8 


0.8497(6) 


0.8224(17) 


0.9678(22) 




8.5985 


0.132427 


12 


0.8441(8) 


0.8198(26) 


0.9712(32) 




8.8323 


0.132169 


16 


0.8410(12) 


0.8203(57) 


0.9754(69) 


1.6950 


7.9741 


0.133179 


8 


0.8411(11) 


0.8079(20) 


0.9606(26) 




8.3218 


0.132756 


12 


0.8340(8) 


0.8154(32) 


0.9777(39) 




8.5479 


0.132485 


16 


0.8328(16) 


0.8076(87) 


0.970(11) 


1.8811 


7.4082 


0.133961 


6 


0.8336(7) 


0.7904(7) 


0.9482(12) 




7.6547 


0.133632 


8 


0.8269(7) 


0.7921(21) 


0.9578(27) 




7.9993 


0.133159 


12 


0.8232(5) 


0.7863(20) 


0.9551(25) 




8.2415 


0.132847 


16 


0.8210(18) 


0.798(12) 


0.972(15) 


2.1000 


7.3632 


0.134088 


8 


0.8172(13) 


0.7778(24) 


0.9518(33) 




7.6985 


0.133599 


12 


0.8097(8) 


0.7757(67) 


0.9579(83) 




7.9560 


0.133229 


16 


0.8091(21) 


0.786(12) 


0.971(14) 


2.4484 


6.7807 


0.134994 


6 


0.8016(9) 


0.7416(17) 


0.9252(23) 




7.0197 


0.134639 


8 


0.7941(9) 


0.7399(37) 


0.9317(48) 




7.2025 


0.134380 


10 


0.7932(9) 


0.7422(37) 


0.9357(48) 




7.3551 


0.134141 


12 


0.7901(13) 


0.7382(54) 


0.9344(70) 




7.6101 


0.133729 


16 


0.7870(9) 


0.756(13) 


0.961(17) 


2.7700 


6.5512 


0.135327 


6 


0.7863(10) 


0.7132(34) 


0.9070(45) 




6.7860 


0.135056 


8 


0.7804(15) 


0.7169(45) 


0.9185(60) 




6.9720 


0.134770 


10 


0.7759(10) 


0.7121(55) 


0.9177(72) 




7.1190 


0.134513 


12 


0.7739(12) 


0.7152(72) 


0.9241(95) 




7.3686 


0.134114 


16 


0.7732(30) 


0.681(34) 


0.880(43) 


3.4800 


6.2204 


0.135470 


6 


0.7573(9) 


0.6558(24) 


0.8659(34) 




6.4527 


0.135543 


8 


0.7501(9) 


0.6560(77) 


0.874(10) 




6.6350 


0.135340 


10 


0.7476(15) 


0.661(11) 


0.885(15) 




6.7750 


0.135121 


12 


0.7430(10) 


0.659(16) 


0.886(21) 




7.0203 


0.134707 


16 


0.7474(13) 


0.642(46) 


0.859(61) 



Table A. 4: Results for the step scaling function 2^ old (in the 'old' scheme). 
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A. 3 Continuum extrapolation of S^ at 

For fixed coupling u the step scaling function defined in Eq. (|A.20f> has a continuum limit, 
cr^ at (u). Neglecting for the moment the uncertainties on the correct values of ct, ct and c^ at , 
we expect the leading-order cutoff effects to be quadratic in the lattice spacing, 

^\u,a/L) = (ri tat (u) + 0(a 2 /L 2 ) , (A.23) 

since 0(a) improvement is employed. Based on this ansatz, Figure [3 in Section |3] illustrates 
the continuum extrapolation of S^ at for a representative subset of our available coupling 
values u = g 2 (L). The coarsest lattices (with L/a = 6) have been omitted from the fits as 
a safeguard against higher order cutoff effects. For the remaining a/L < 1/8, the one-loop 
cutoff effects are quite small, see Figure E3 

Although these extrapolations are entirely compatible with an approach to the contin- 
uum limit quadratic in a/L, we also have investigated extrapolations linear in a/L. These 
as well yield reasonable fits with consistent results and even comparable total x 2 /dof (when 
summing up the % 2 -s belonging to the individual fits at the ten it-values) so that the form 
of the lattice spacing dependence can not be decided on the basis of the data. Therefore, we 
have studied the influence of the imperfect (i.e. only perturbative) knowledge of some of the 
improvement coefficients in more detail. 

Since the usage of the two-loop approximation (jA.21j) for ct in the calculation of the 
correlation functions (and thereby also in the step scaling function S^ at ) should cancel the 
main contributions from the related boundary terms, we only address its effect originating 
from the fixing of the renormalized coupling, the values of which were taken over from 
Ref. ^3] with c t still set to one- loop. Changing c t from one- to two- loop also in this step 
then requires to adjust the bare coupling and the value of the critical quark mass accordingly 
before the simulations for S^ at can be performed. We have done this analysis for the largest 
fixed coupling, u = 3.48, where the uncertainty in ct is largest and thus its effect most 
pronounced. At L/a = 6 the resulting change in E^ at turns out to lie clearly inside the 
statistical errors, and this effect will even get smaller for decreasing a/L. On the other hand, 
if we just compute S^ at with the one-loop value of ct as in in the computation of the coupling, 
we found the results, now for u = 2.77 and L/a = 6, 8, to be indistinguishable within errors, 
too. We conclude that any small uncertainty present in ct beyond the available two-loop 
estimate is numerically unimportant for the cutoff dependence of S^ at . 

Regarding the 0(a) improvement coefficient ci, we followed the same line as in ^3] and 
assessed its influence on our results by artificially replacing the one-loop coefficient in the 
expression (|A.22|) by ten times its value. I.e. we set ci to £%' = 1 — 0.180 g$ in some additional 
simulations at u = 3.48, and the outcome is that the corresponding estimates on S^ at for 
L/a = 6 still differ by around 1.5%, while for L/a = 8 they already agree within their 
statistical errors. As this difference drops further for growing L/a and/or smaller couplings, 
a possible imperfection of c t does not affect the results on £ J at either. 

Finally, we also checked for the influence of the 0(a) improvement coefficient in the 
static-light axial current, c^ at , by analyzing our data with = instead of the one-loop 
value (|A.9|) . Whereas the related change in Z| tat is of the order of a few percent and hence 
still substantial, it largely cancels in the ratio of Eq. (|A.20|) so that this effect is no more 
significant for given its statistical errors. 

All in all these findings demonstrate that at the level of our precision linear a-effects 
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in the data on E^ at are negligible, and extrapolations using (a/L) 2 -terms as the dominant 
scaling violation are justified indeed. 

Appendix B Perturbation theory 

This appendix provides a few details on the perturbative computations, which were required 
to obtain the one-loop expression for h(d/L, u), Eq. (|2.12|) . the two-loop anomalous dimension 
and the one-loop estimates of the discretization errors of the step scaling function S^ at . Note 
that here we restrict ourselves to the case of the modified (or 'new') scheme introduced via 
Eq. (|2~T5|) in this paper, because the perturbation theory of the original scheme defined 
through Eq. (|2.7j) has been extensively discussed already in Ref. JJj where also more details 
on the different steps involved can be found. 

The correlation functions /f tat , /Ia* and f\ are expanded in powers of the coupling 
as explained in J2j and |H2] , and the analogous expansion of /j 1 * 1 is explained below. 



B.l The correlation function / hh 



i 

hh 



After integrating out the static quark fields, the correlation function /j 1 can be written as 
/1V3) = |2 Yl (te{u(x,0)U(x + aQ,0)---U(x+(L-a)6,0) 

XI, X2 

xU((L-a)6,0)- 1 U((L-2a)6,0)~ 1 ---U(0,0)- 1 } V (B.l) 

J xn=0 / 



where the trace is taken over colour indices only. 

Writing U(x,/jL) = exp {g^aq^x)}, with the gluon field q^x) = q®(x)T a , where T a are 
the antihermitian generators of the gauge group, the function /j 111 can be expanded in the 
bare coupling, 

ft\d) = 3 + g 2 fi Hl) (d) + O(g 4 ), d = \x 3 \ . (B.2) 

Determining the one-loop coefficient J^h.C 1 ) am0 unts to calculating and summing the 
diagrams shown in Figure IrTTI using the gluon propagator D^ u (x,y) (with SF boundary 
conditions) given in Ref. 





Figure B.l: One-loop diagrams contributing to /j 1 . The two dots on the left are at xq = 0, 
the dots on the right are at xq = L. 
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For comparison with the non-perturbative results shown in Figure 0J Subsection 12.21 we 
consider the one- loop coefficient h^ l \d/L) of Eq. (j2.13|) . which reads 

ffi-Hd/L) = \ {fi h ' (1 \d) - / 1 hh ' (1) (L/2)} , (B.3) 

and which we can obtain analytically. (For the other quantities considered in this appendix, 
the diagrams are calculated numerically.) Only the last of the diagrams in Figure IB. II 
contributes, and we thus can write 

h {1) (d/L) = ^rj Yl X) { D oo(x,y)\ X3 =d, 3/3=0 ~ Doo(x,y)\ X3=L/2i ya=0 } 

X0,V0 XltVl %2,V2 

= ^3^2{ eiP3d - emL/2 } E^fo'^CO.O.j*)), (B.4) 
P3 £0,2/0 

with the momentum-space gluon propagator doo(xo,yo;p) defined in Ref. jSH]. The p$ = 
term does not contribute to the sum, and using the explicit form of doo(xo, yo; p) for p ^ 0, 
one can show that 

a 2 ^ d 00 (x ,y ; (0,0,p 3 )) = (B.5) 
%o,yo ^ 3 

with f>z = \ sm (^) • Thus we see that h^ l > (d/L) is just given in terms of the one-dimensional 
scalar propagator on a periodic lattice with length L, which has no lattice artifacts. This 
eventually leads to the result quoted in Eq. Q2.14JI . The absence of any lattice spacing 
dependence is a consequence of the special kinematics, namely the summation over xi,X2, 
but will of course not be exact in higher orders of perturbation theory. 



B.2 Anomalous dimensions 

In order to precisely connect to the RGI current, it is important to obtain the anomalous 
dimension of the static-light axial current in the SF scheme at two-loop order. The anomalous 
dimension is expanded as 

i(g) = -g 2 { 70 + 7i SF 9 2 + (V) } , (B.6) 

with 70 = — l/(47r 2 ) and the two- loop anomalous dimension in the SF scheme, 7^ . 

With the perturbative expansions for the various correlation functions, the ratio Hi of 
Eq. (|A.16|) can be written as a series 

3i(go, a/L) = Z[°\a/L) + g 2 ^{a/L) + O(g A ) . (B.7) 

Accordingly, this allows us to expand the SF renormalization constant Z^ at = St /Si (see 
Eq. (Oty ) as 

^stat = z stat,(0) + 5 2^tat,(l) + Q ^ ( g g) 

With the one-loop relation between the bare lattice current and the renormalized static 
axial current in the MS scheme, the anomalous dimension in the MS scheme, Eq. Q5.5[) 38 , 
I39| l4()j. can be converted into the SF scheme. The renormalization constant relating the SF 
scheme and the MS scheme is obtained from the relation between the SF scheme and the 
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bare lattice current, Eq. ()B.8|) . the connection of the bare lattice current and a 'matching 
scheme' [HUEl] and the relation between the latter and the MS scheme [SOI- Here the 
matching scheme is defined by the requirement that the renormalized static-light axial current 
at scale \x = mh equals the relativistic axial current with a heavy quark mass m\ up to terms 
of 0(l/mh), and the current in the relativistic theory is normalized by current algebra 
(imposing the chiral Ward identities). 6 Following the steps in Ref. J21> this analysis finally 
yields 

(9 = 0.0: 7! SF = — i-rr { 0.22(2) - 0.0552(13)iV f } , (B.9) 

= 0.5: = t—^ {0.10(2) -0.0477(13)iVf}, (B.10) 

(47r) z 

e = 1.0 : 7? F = t—^ { -0.08(2) - 0.0365(13)iV f } . (B.ll) 
B.3 Discretization errors 

The one-loop expansion at hand is also helpful to study discretization errors in the step 
scaling function. Using Eq. ()B.8|) . the step scaling function at lattice spacing a is expanded 
as 

Sj at (n, a/L) = 1 + u £^ at ' (1) ( a /L) + 0(u 2 ) , (B.12) 
and its continuum limit cr^ at (n) as 

stat/ \ 1 , stat,(l) 2 Stat, (2) 3 A /T3 -. „n 

a A [u) = l + ua A y '+ua A + 0(u ) , (B.13) 

with 

stat,(l) , /„x 

<7 A = ln(2) 70 , 

ff rtat,(2) = l ln 2 (2)7 2 + ln 2 (2)6o7o + ln(2)7i> (B.14) 

As a measure for the discretization errors, we define 

(B , 5) 

with a perturbative expansion 

S(u,a/L) = 5 {1) (a/L)u + 0(u 2 ). (B.16) 

The one-loop coefficient 5^ versus the lattice spacing squared is shown in Figure IB.2I for 
different values of 9. In the range of lattice spacings where our non-perturbative calculation 
is performed, the discretization errors at one-loop level are smaller than 1% x u, giving rise 
to the hope that also the non-perturbative discretization errors are reasonably small. A 
welcome feature of 8^' at 8 = 0.5 is that it is entirely dominated by the leading a 2 /L 2 -term 
in the a-expansion. 

6 Since the wording in Ref. |17| is not completely clear on this, we point out that in that reference 

refers to what we call the MS scheme here as well as in Ref. |17|. 
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Figure B.2: Discretization errors of the step scaling function at one-loop level. The con- 
tinuum extrapolation of the non-perturbative results uses (a/L) 2 < 0.02. 



Appendix C Continuum step scaling function and matching at L = 1.436 r 

In this last appendix we briefly discuss our parametrization (i.e. the interpolating fit) of the 
continuum step scaling function and some details on the calculation of Z| tat at the matching 
scale 2L max = 1.436 r^. 

C.l Fits and error determination in the scale evolution 

As described in Appendix lA.3| the continuum step scaling function cr^ at (-u) has been obtained 
by extrapolating the lattice data on S^ at (u,a/L) at fixed u to the continuum limit. The 
next step is now to solve the recursion specified through Eqs. ()3.1j) - l|3.4|) . In practice this 
is done by first representing the results for cr^ at (tt) in Tabled by a fit and then solving the 
recursion, Eq. 1)3.4)) . with (u) given by the fit function. 

Guided by the analysis for the step scaling function of the pseudoscalar density crp in 
Ref. ^3] an d the perturbative expansion discussed in Appendix IB. 31 

af^(u) = 1 + s u + s x u 2 + s 2 u 3 + ... + s n u n+1 (C.l) 

is chosen as fit ansatz. The two non-trivial leading terms are restricted by perturbation 
theory, 

stat.(l) stat,(2) /.-. 

so = cr A , si = a A (C.2) 

cf. Eqs. (|B.14)) . Up to three additional free fit parameters were allowed for. All of these 
fits represent the data in Table Q well, and we decided to quote the two-parameter fit (the 
curve of which is shown in Figure EJ) as the final result for the functional form of <r^ at . To 
check that the polynomial fits are stable, we also investigated fits where only sq or even no 
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coefficient at all is constrained to its perturbative value. This leads to consistent results for 
a^ at («); particularly the latter fit then reproduces the perturbative prediction for so- 

Having chosen a definite expression for cr^ at (u) , the solution of the associated recursion 
is unique. Since the errors on the step scaling function stem from different simulation runs 
and are hence uncorrelated, the errors on the fit parameters in the polynomial (|C.1|) for 
cr^ at (u) and those on the Vk-s calculated from it can be estimated straightforwardly by the 
standard error propagation rules. Finally, by increasing the number of free fit parameters 
(while fixing so, si to perturbation theory) as mentioned above, we convinced ourselves that 
the systematic error induced by the choice of fit functions is well under control: in fact, we 
then observed the expected pattern of finding slightly different errors but compatible results 
at comparably good overall fit quality. 

C.2 Calculation of Z^ tat at the low-energy matching scale 

The total renormalization factor Zrgi introduced in Section |1] involves the value of the 
renormalization constant Z| tat at our particular matching point: Z^ a,t (go, L/a) | i=1 436 . 
As the latter connects a bare matrix element of the static-light axial current to the one 
renormalized in the SF scheme, this amounts to calculate Z^ at for a range of bare couplings 
commonly used in simulations in physically large volumes. 

To extract Z^ at we exploit the fact that the required pairs (L/a, j3) that match the 
condition L/a = 1.436 r§/a have already been determined for the relevant /3-range in Ap- 
pendix C of 13 by utilizing the known parametrization of ln(a/ro) in terms of f3 from 
Ref. |34| . We thus could take over these pairs and computed Z^ tat for 9 = 0.5 from the 
renormalization condition ()A.19|) at the corresponding values K = n c of the critical hopping 
parameter jHS]- The results for Zj£ (go,L/a)\ L=li36ro using the one- and two-loop expres- 
sions for ct, cf. Eq. (|A.21|) . are given in Table RTT1 The difference originating from the two 
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Z A 
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6.0219 


0.13508 


0.6926(15) 


0.13504 


0.6932(10) 


10 


6.1628 


0.13565 


0.6810(17) 


0.13564 


0.6824(16) 


12 


6.2885 


0.13575 


0.6786(18) 


0.13574 


0.6795(16) 


16 


6.4956 


0.13559 


0.6777(17) 


0.13558 


0.6763(18) 



Table C.l: Results for Z^ 8 * (go, L/a) at fixed scale L = 2L max = 1.436 ro with c t being 
set to either one- or two-loop. The critical n- values k = k c \65\ are the same as used for 
Table C.l of Ref . JT^. 

perturbative approximations for ct is completely covered by the statistical errors so that we 
again consider c\ ~ loop to already account for the dominant part of the boundary cutoff effects 
in the gauge sector. Similarly to the discussion in Appendix IA. 31 the influence of the bound- 
ary improvement coefficient c t in the fermionic sector can also be neglected at the level of our 
precision. The parametrization of the results for Z^ at (go, L/a)\ L=1 43 g ro by a polynomial fit 
in ((3 — 6), with c^~ loop and c^ at from one-loop perturbation theory, is quoted in Sectional 
where the coefficients in the first block of Tableware to be combined with Eq. (|4.1|) . The 
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smooth dependence of Z s ^ &t on (3 in the studied region of bare couplings suggests that this 
representation can also be slightly extended down to (3 = 6.0 (even though we could not 
directly simulate that point for the same reason as in case of Zp |13|). and we therefore 
regard it as a reliable representation of our data over the whole range 6.0 < (3 < 6.5. 

As in Appendix lA.31 we also set c^ at = instead of one-loop in the analysis of the data on 
■^A tat (5o ; L/a)\ L=l 436r , and the results are listed in the middle part of Table RT2l In contrast 



L/a (3 



6/5o 2 



»7Stat 



K 



4 


5.6791 




6 


5.8636 




8 


6.0219 


0. 


10 


6.1628 


0. 


12 


6.2885 


0. 


16 


6.4956 


0. 



13508 0.6736(14) 

13565 0.6633(17) 

13575 0.6621(17) 

13559 0.6627(17) 



0.15268 0.6923(13) 

0.15451 0.6315(16) 

0.15341 0.6075(13) 

0.15202 0.5964(18) 

0.15078 0.5971(32) 

0.14887 0.5991(52) 



Table C.2: Results for Zf at (g ,L/a) at fixed scale L = 2L max = 1.436 r for cf at = and 
unimproved Wilson fermions (i.e. c sw = and thus, ca = = too), where c t was kept 
at its two-loop value. 

to the step scaling function, which did not change appreciably under this replacement of the 
value for c A iat , we observe an effect of about 3% in Z^ at at the low-energy matching scale 
L = 2L max = 1.436 r . 

Furthermore, we addressed the case of unimproved Wilson fermions by also setting 
c S w = in the relativistic fermion action and, after having computed the needed estimates of 
the critical hopping parameter for this situation, carried out the additional runs to determine 
the renormalization constant. In this case the pairs {L/a, [3) were extended to lower values 
of (3 in order to be able to make contact with the /3-region that is typically employed in 
simulations to calculate the bare matrix element defining the B-meson decay constant, as 
e.g. those in Refs. jlOUllj . The resulting estimates on Z^ at are shown in the right part of 
Table EZ3 and the corresponding polynomial representations for both aforementioned cases 
are as well found via Eq. (|4.1j) together with the two lower blocks of Table |3 

We conclude this discussion with the general remark that the uncertainties of the en- 
tering critical K-values (of 1-2 and 2-4 on the last decimal place in the case of c sw = 
non-perturbative and c sw = 0, respectively) do not affect the Z-factors significantly. 
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